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We suggest a new Hamiltonian lattice approach, using a regularisation motivated by deep inelastic scattering. 
We discuss the relation between distribution functions and the f\ structure function. We have tested the method 
by computing the critical behaviour of the scalar model and find agreement with scaling behaviour and with 
results by Liischer and Weisz. 



The computation of quark distribution func- 
tions, structure functions and of the hadronic 
mass spectrum directly from QCD requires non- 
perturbative methods. The most important non 
perturbative approach to compute these quanti- 
ties is lattice gauge theory. Lattice gauge the- 
ory can be formulated in both Euclidean space 
or Minkowski space. Here we shall explore the 
usefulness of the Hamiltonian or transfer-matrix 
approach in Minkowski space as an alternative 
to the commonly used statistical methods in Eu- 
clidean space 1 . From a numerical point of view, 
however, it is extremely difficult to actually di- 
agonalise a Hamiltonian ( or transfer matrix) 
even on very small lattices with the currently 
available numerical techniques. Euclidean lat- 
tice gauge theory is up to now the only non- 
perturbative method which is capable of nu- 
merically describing the low-energy part of the 
observed hadronic mass spectrum directly from 
QCD (Glueball masses, however, have already 
been computed in the Hamiltonian approach[l]). 
But unfortunately, the computation of distribu- 
tion functions or structure functions is very diffi- 
cult in Euclidean space. So far, the first moments 
of valence quark structure functions can be mea- 
sured. Important work in this direction has been 
done by Schierholz and co-workers [2]. The com- 
putation of the first moments of the sea quark 
contribution to the structure functions is feasible 
in principle but in practice very complicated and 

1 Both approaches should be equivalent if the measured 
correlation length is much greater than the lattice spacing 
For a small correlation length, a "renormalisation of the 
speed of light" can be necessary. 



has not been done yet. A computation of quark 
distribution functions beyond the first moments is 
still beyond reach of the currently available meth- 
ods. The computation of distribution functions 
(not only the first moments thereof) would be 
relatively simple if it were possible to obtain the 
hadronic wave functions from a Minkowski space 
calculation. Two other problems of the Euclidean 
approach 

- finite density thermodynamics and 

- excited states of mass spectra 

are expected to be simpler in the Minkowskian 
approach, if it were possible to diagonalise the 
transfer matrix of the field theory under inves- 
tigation. A further advantage of this approach 
is that it is closer to the intuitive picture of a 
Hadron as a many body bound state of quarks 
and gluons with wave functions describing the 
probability amplitude of finding a certain config- 
uration of quarks and gluons inside the proton. 
It is, however, this very advantage which causes 
the numerical problems because a hadron consists 
of a large avarage number of constituents which 
diverges in the thermodynamic limit. A many 
particle system is very difficult to treat numeri- 
cally and even more so a system with fluctuating 
(non conserved) constituent number. To render 
things worse, mean field many body techniques 
(e.g. Hartree-Fock) are not viable for a relativis- 
tic field theory in fewer than 4+1 space-time 
dimensions. 



So how could the QCD Hamiltonian be diag- 
onalised with reasonable numerical effort? Our 
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idea is to consider experiment as a guide-line. 
More specificly we consider deep inelastic scat- 
tering experiments where the composite nature 
of hadrons becomes visible. Our aim is to de- 
scribe only properties of the constituents in the 
the hadron than we are actually able to measure. 
This has three immediate consequences. Firstly, 
the observed quarks depend on the experimen- 
tal resolution, according to the renormalisation 
group improved parton model. They are effective 
particles, not elementary ones. Quarks observed 
at low resolution consist of quarks and gluons at a 
higher resolution scale. Hence the lattice spacing 
has to be chosen to be of the same order of mag- 
nitude as the experimental resolution. Secondly, 
the resolution of the experiment (= the wave- 
length = l/\q] of the exchanged virtual photon 
with four-momentum q ) is not Lorentz invariant 
and so are the effective quarks. In the Breit-frame 
(where the hadron four-momentum takes on the 
form P = (E, 0.0, P3) and where the momentum 
of the virtual photon q = (0,0,0,— Q)) the in- 
verse resolution coincides with the relativistic in- 
variant Q := \/—q 2 - In this frame, the Bjorken 
scaling variable x := can be interpreted as 

the fraction x = of the momentum p of the 
struck parton inside the total proton momentum 
P3. Therefore we chose this frame as the lattice 
rest-frame. Finally and most importantly the ob- 
servable parton momenta necessarily lie in the in- 
terval < P3 < P3 if the hadron is stable. Hence, 
for every finite lattice regularisation the number 
of observable states in the Fock space is finite. 
The mathematical reason for this is that only a 
finite number of particles with discrete momenta 
can share a given total momentum. The infinite 
momentum frame shares a similar feature. There 
is, however, a major advantage of using a high 
momentum frame which is not the infinite mo- 
mentum frame: We have a finite space — like lat- 
tice spacing at our disposition which is not possi- 
ble in the front-form nor is it possible with infinite 
hadron momentum(!) This is an important fea- 
ture because without a finite lattice spacing, no 
dimensionless correlation length £ = — - — can be 
defined from the hadron mass uir. In the con- 
tinuum limit £ — y 00, however, we end up with 



an infinite hadron momentum since the hadron 
velocity v = P 3 /E = Psa/^/l/C 2 + P^a 2 -> 1 
approaches the velocity of light. The same thing 
holds for the Bjorken limit Q — » 00, x = const, 
since in the Breit frame the parton momenta 
p 3 = Q/2 — > 00 as well as the proton momen- 
tum P3 = Q/2/x — > 00, diverge in this limit. 

As a first test of the Breit method, we have 
applied it to the (j> 4 model in four space-time di- 
mensions. The Lagrangian density of this model 
is 

In Euclidean lattice calculations, the mass m and 
the coupling constant g are usually replaced by 
the hopping parameter k and the lattice coupling 
constant A as follows m 2 = 1-2A —8,0 = 
This model shows two phases: the asymmetric 
phase and the symmetric one. We considered the 
latter. As one approaches the critical line in the 
k,\ parameter space, the correlation length £ = 
— - — diverges, where mo is the mass of the ground 
state. Close to the critical line, the ground state 
masses should obey the following universal scaling 
law: 

l/£ = m R a ~ CY2|logr|- 1/6 

where r = K ~ K|1 and k c is the critical couplong 
where the correlation length diverges. Fig.[l] 
compares our results to semianalytic calcula- 
tions[3] showing that the asymptotic scaling law 
is nicely reproduced. In Fig. [2] the mass ratios of 
the spectrum are displayed. We observe scaling 
close to the critical point. 

One comment on the scalar (j> 4 theory is in or- 
der: If we had used the zero-momentum sector 
where the bound state is at rest, the problem 
would have been a complicated many-body prob- 
lem. The 'parton cloud' would have appeared [4]. 
In a fast-moving frame, however, the one-particle 
sector dominates the bound-state. 

We have designed our method in view of QCD. 
Therefore, it is useful to calculate the structure 
functions in our approximation. A straightfor- 
ward calculation using the impulse approximation 
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yields the following expression 
F 1 (x,Q) = J2el/Hf { ; ) (x,l 



^H/i'W,^)} 



where fs (p3, P3, 



J d 2 PL B(p)wi i} (p,P 3 ,fi) 

and where wy(p,P 3 ,n) := < ps\ps> 

is the probability of finding a quark with helicity 
s and momentum p inside a hadron with four- 
momentum P and Pauli-Lubanski spin vector S. 
In analogy, w and / refers to an antiquark. The 
quantity B is defined as B(p) := -5- / 

2 p 3 \/m 2 +p 2 

If we neglect parton masses and transverse mo- 
menta in the Bjorken limit we end up with the 
well known interpretation of the F\ structure 
function as a linear combination of quark distri- 
bution functions q s ,i(%,Q) = fs l \x, l,2x fi/Q) = 
Psfs^ (p3, P3, A*)- The same holds for F2 and g\. 



Fig.l The ground state mass tjir in lattice units 
(a = 1) versus k for A = 0.00345739 (A = 0.01 
in Ref.[3]). The dots correspond to results of 
Ref.[3]. Our results correspond to a 7 3 lattice 
(dashed line) and to a 9 3 lattice (solid line). 



In conclusion, the diagonalisation of a relativis- 
tic Hamiltonian is much simpler in a large mo- 
mentum sector than in the zero momentum sec- 
tor. We have tested this for the case of the (f) 4 - 
model where we obtained asymptotic scaling in 
good agreement with the semianalytical results of 
Liischer and Weisz. In deep inelastic scattering 
this corresponds to taking into account partons 
with experimentally observable momenta only. 
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Fig. 2 The lowest lying mass spectrum versus k. 
The ground state mass is set to one. A as in 
Fig.[l]. 
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